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Abstract 

It is well known that a supercritical single- type Bienyame-Galton- 
Watson process can be viewed as a decomposable branching process 
formed by two subtypes of particles: those having infinite line of de- 
scent and those who have finite number of descendants. In this paper 
we analyze such a decomposition for the linear-fractional Bienyame- 
Galton- Watson processes with countably many types. 
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1 Introduction 

The Bienayme-Galton- Watson (BGW-) process is a basic model for the stochas- 
tic dynamics of the size of a population formed by independently reproducing 
particles. It has a long history [2] with its origin dating back to 1837. This 
paper is devoted to the BGW-processes with countably many types. One 
of the founders of the theory of multi-type branching processes is B.A. Sev- 
astyanov [8], [9]. 
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A single-type BGW-process is a Markov chain {Z^}^ =0 with countably 
many states {0,1,2,...}. The evolution of the process is described by a 
probability generating function 

oo 

f(s) = Y,PkS k , pi<l, (1) 

where pk stands for the probability that a single particle produces exactly k 
offspring. If particles reproduce independently with the same reproduction 
law (jl|, then the chain {Z^ n '}^ =0 represents consecutive generation sizes. In 
this paper, if not specified otherwise, we assume that = 1, the branching 
process stems from a single particle. Due to the reproductive independence 
it follows that f^(s) = E(s z<n> ) is the n-th iteration of f(s). 

Since zero is an absorbing state of the BGW-process, = F(Z^ = 0) 
monotonely increases to a limit q called the extinction probability. The latter 
is implicitly determined as a minimal non-negative solution of the equation 

f(x) = x. (2) 

A key characteristic of the BGW-process is the mean offspring number M = 
f'(l). In the subcritical (M < 1) and critical (M = 1) cases the process is 
bound to go extinct q — 1, while in the supercritical case (M > 1) we have 
q < 1. Clearly q = if and only if p = 0. 

In the supercritical case the number of descendants of the progenitor 
particle is either finite with probability q or infinite with probability 1 — q. 
Recognizing that the same is true for any particle appearing in the BGW- 
process we can distinguish between skeleton particles having an infinite line 
of descent j6] and doomed particles having a finite line of descent. Graphically 
we get a picture of the genealogical tree similar to that given in Figure 1. 

If we disregard the doomed particles, the skeleton particles form a BGW- 
process with a transformed reproduction law excluding extinction 

/w = /(s(1 7^ 9) ' 9 (3) 

and having the same mean M — M > 1. Formula ^ is usually called the 
Harris-Sevastyanov transformation. On the other hand, the doomed particles 
form another branching process corresponding to the supercritical branch- 
ing process conditioned on extinction. The doomed particles produce only 
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Figure 1: An example of a BGW-tree up to level n = 10. Solid lines repre- 
sent the infinite lines of descent and dotted lines represent the finite lines of 
descent. 

doomed particles according to another transformation of the reproduction 
law f(s) = f(sq)/q, which is usually called the dual reproduction law and 
has mean M = f'(q) < 1. The supercritical BGW-process as a whole can be 
viewed as a decomposable branching process with two subtypes of particles 
[Tj Ch. 1.12]. Each skeleton particle must produce at least one new skeleton 
particle and also can give rise to a number of doomed particles. In Section 
[2] we describe in detail this decomposition for the single type supercritical 
BGW-processes. 

In the special case when the reproduction generating function ([I]) is linear- 
fractional many characteristics of the BGW-process can be computed in an 
explicit form [5]. In Section [3] we summarize explicit results concerning de- 
composition of a supercritical single-type BGW-processes. 

Section [4] presents the BGW-processes with countably many types. Our 
focus is on the linear-fractional case recently studied in [7] . The main results 
of this paper are collected in Section [5] and their derivation is given in Section 
|6j The remarkable fact that a supercritical branching process conditioned on 
extinction is again a branching process was recently established in jl] in 
a very general setting. In general, the transformed reproduction laws are 
characterized in an implicit way and are difficult to analyse. This paper 
presents a case where the properties of the skeleton and doomed particles 
are very transparent. 
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2 Decomposition of a supercritical single-type 
BGW-process 

The BGW-process is a time homogeneous Markov chain with transition prob- 
abilities satisfying 

oo 

£^ = (/ (n W- 

In the supercritical case with mean M > 1 and extinction probability q < 1, 
using the property Y^JLo PijQ* — Q\ we can get another set of transition 
probabilities putting 



p(«) ._ p( n ) n j-i 

The transformed transition probabilities also possess the branching property 



oo 



j=0 

where f^ n \s) = f ( - n \sq)/q is the n-th iteration of the so-called dual generat- 
ing function 

f(s) = = f>s fe , p k =p k q k -\ k > 0. 

The corresponding dual BGW-process is a subcritical branching process 
with offspring mean M = f'(q) < 1, see Figure [ij The dual BGW-process is 
distributed as the original supercritical BGW-process conditioned on extinc- 
tion: 

p(zw = j\z<» = izw = o) = nz^=j,z^ = o\z^ = ^ 

[ 31 ' UJ P(ZM = 0|Z(°) = z) 

= g-T(Z (oo) = 0|Z (n) = j)i^ n) 



p( n ) 

ij " 
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Figure 2: Duality between the subcritical and supercritical cases. Left: a 
supercritical generating function ([I]) with two positive roots (q, 1) for the 
equation Right: the dual generating function drawn on a different scale. 

The two parts of the curve on the left panel of Figure [2] represent two 
transformations of the supercritical branching process. The lower-left part 
of the curve, replicated on the right panel of Figure [2] using a different scale, 
gives the generating function of the dual process. The upper-right of the 
curve on the left panel corresponds to the Harris- Sevastyanov transforma- 
tion The function ^ is the generating function for the probability 
distribution 

Po = 0, p fe = f>/* V*(l - 

i=k ^ 

with the same mean M = M as the original offspring distribution. It is easy 
to see that the n-th iteration of f(s) is given by 

Kn) ( x _ f M (s(l-q) + q)-q 
1 [) ~ l-q 

Looking into the future of the system of reproducing particles we can 
distinguish between two subtypes of particles: 

• skeleton particles with infinite line of descent (building the 
skeleton of the genealogical tree), 

• doomed particles having finite line of descent. 

These two subtypes form a decomposable two-type BGW-process {(S^ n \ D^)} 
with + = Z^ n \ The joint reproduction law for the skeleton particles 
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has the following generating function 

F(s,t) = E(s*V (1) ) = /(s(1 - g) i + t J ) ~ /(tg) . 
A check on the branching property for the decomposed process is given by 

TRY Z^J-Z^ ■] \ 

E(s s(n) t D(n) \z^ > o) = {s 1 / ; h^ll 

V I > p(^(oo) > 0) 

_E(s s(n) t D(n) )-E(t D(n) l {z ^ =0} ) 



l-q 

f in) (s(l-q)+tq)-fW(tq) 
l-q 



= F^ n \s,t). 



The original offspring distribution can be recovered as a mixture of the joint 
reproduction laws of the two subtypes 

/W( S ) = (l-g)FW( S , S )+g/»( S ). 

Observe also that the total number of offspring for a skeleton particle has a 
distribution given by 



f(s) = F(s, s) = /(S) 1 _ f(Sq) = f>A p k = (1 - q)-\l - q k ) Pk , 

Q k=0 

with mean M = It follows, 

l-q ' 



M M ■ (M M) | - 



Q 

Q 

and we can summarize the relationship among different offspring means as 

M <1 < M = M < M. 

3 Linear-fractional single-type BGW-process 

An important example of BGW-processes is the linear-fractional branching 
process. Its reproduction law has a linear-fractional generating fuction 

M = ho + — — (4) 

1 + m — ms 
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fully characterized by two parameters: the probability h = p of having no 
offspring, and the mean m of the geometric number of offspring beyond the 
first one. Here h\ — 1 — ho stands for the probability of having at least one 
offspring. Notice that with h = — L, the generating function Q describes a 
Geometric (j^) distribution with mean m. If ho = the generating function 
^ gives a Shifted Geometric (jt^) distribution with mean m + 1. If m = 0, 
we arrive at a Bernoulli (hi) distribution. 

Since the iterations of the linear-fractional function are again linear- 
fractional, many key characteristics of the linear-fractional BGW-processes 
can be computed explicitly in terms of the parameters (ho, m). For example, 
we have M = hi(l + m), and if M > 1, we get 

, , 1 + m- M 

q = h (l + m ) = . 

m 

The dual reproduction law for ^ is again linear-fractional 

f(s) = h + - hlS ; , ho= m - , m = h /h 1 , 
1+m — ms m + 1 

with M = 1/M. The Harris-Sevastyanov transformation in the linear- 
fractional case corresponds to a shifted geometric distribution 

~ s 

f(s) = — , m = m(l — q) = M — 1. 

1 + m — ms 

Interestingly, the joint reproduction law of skeleton particles 
hi f s(l — q) +tq tq 



F(s,t) 



1 — q \l + m — m(s(l — q) +tq) 1 + m — mtq 
s 1 



1 + m — m(s(l — q) + tq) 1 + m — mt 
1 1 



1 + m — ms — (m — m)t 1 + m — mt 
has three independent components: 

• one particle of type 1 (the infinite lineage), 

• a Geometric (jt^) number of offspring each choosing inde- 
pendently between the skeleton and doomed subtypes with 
probabilities m/m and (m — m)/m, 
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• a Geometric (j^) number of doomed offspring. 

Observe that even though both marginal distributions f(s) and f(s) are 
linear-fractional, the decomposable BGW-process (S^ n \D^) is not a two- 
type linear-fractional BGW-process. The distribution of the total number of 
offspring for the skeleton particles in not linear-fractional 

m = — — 1 ; 

1 + m — ms 1 + m — ms 

and has mean 

M = l+ m + m = M+(M + l)m. 



4 BGW-processes with countably many types 

A BGW-process with countably many types 

ZW = (ZW^...), n = 0,1,2,... 

describes demographic changes in a population of particles with different 
reproduction laws depending on the type of a particle % G {1,2, . . .}. Here 
z\ n ^ is the number of particles of type % existing at generation n. In the 
multi-type setting we use the following vector notation: 

x = (xi, x 2 ...), 1 = (1, 1 . . .), e i = (l { i=i } , l{i= 2 >, ■■■), 
xy = (x 1 y 1 ,x 2 y2,---), x -1 = (xi 1 ,a;J 1 ,...), x y = xf • • • , 

we write x*, if we need a column version of a vector x. 

A particle of type i may produce random numbers of particles of different 
types so that the corresponding joint reproduction laws are given by the 
multivariate generating functions 

/,(s)=E(s z(1) |Z(°)=e,). (5) 

The offspring means 

M lJ =E(zj 1) |Z(°) = e l ) 

are convenient to summarize in a matrix form M = (Mjj)°° =r For the n-th 
generation the vector of generating functions f ^ (s) with components 

^(s)=E(s z(n) |Z(°) = e i ) 
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are obtained as iterations of f(s) with components ([5]), and the matrix of 
means is given by M n . The vector of extinction probabilities q = (51, ?2, • • •) 
has its i-th component qi defined as the probability of extinction given that 
the BGW-process starts from a particle of type %. The vector q is found as 
the minimal solution with non-negative components of equation f(x) = x, 
which is a multidimensional version of ([2]). 

From now on we restrict our attention to the positive recurrent (with re- 
spect to the type space) case when there exists a Perron-Frobenius eigenvalue 
p for M with positive eigenvectors u and v such that 

vM = pv, Mu' = pu\ vu* = vl* = 1, p- n M n ->■ uV, n -> oo. 

In the supercritical case, p > 1, all < 1 and we can speak about the 
decomposition of a supercritical BGW-process with countably many types: 
(S^D^™)). Now each type is decomposed in two subtypes: either with 
infinite or finite line of descent. The decomposed supercritical BGW-process 
is again a BGW-process with countably many types whose reproduction law 
is given by the expressions 

F|(g|t)= /.(■(! -q)+tq)-/«(tq) / <(t) = *M. 

Linear-fractional BGW-processes with countably many types were stud- 
ied recently in [7]. In this case the joint probability generating functions ^ 
have a restricted linear-fractional form 

fi(s) = h i0 + — 3 ~ l JiJ . (6) 

i + m-m2^ j=1 gjSj 

The defining parameters of this branching process form a triplet (H, g, m), 
where H = (^)^ =1 is a sub-stochastic matrix, g = (</i, <?2> • • •) is a proper 
probability distribution, and m is a positive constant. The free term in ^ 
is defined as 

oo 

hio — 1 — hy. 

3=1 

As in the case of finitely many types jl], the denominators in ^ are neces- 
sarily independent of the mother type to ensure that the iterations are also 
linear-fractional. This is a major restriction of the multitype linear-fractional 
BGW-process excluding for example decomposable branching processes. 
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It is shown in [7] that in the linear-fractional case the Perron-Frobenius 
eigenvalue p, if exists, is the unique positive solution of the equation 



OG 

m^p^gH^l 1 = 1. (7) 

k=l 

In the positive recurrent case, when the next sum is finite 

oo 

(3 = m J2kp- k gH k l\ (8) 



k=l 



the Perron-Frobenius eigenvectors (v, u) can be normalized in such a way 
that vu 1 = vl* = 1. They are computed as 

oo 

u* = (1 +m)/r 1 ^p- fe H fe l t , (9) 



k=l 



v = Tz^5>~ fe s Hfe - ( 10 ) 

1 + m ^— ' 



" " k=0 

In the supercritical positive recurrent case with p > 1 and < oo the 
extinction probabilities are given by 

q= 1 - (p- 1)(1 +m)- 1 (3u. (11) 

Observe that gu* = and 

gq* = ^. 12 

m 

The total offspring number for a type i particle has mean 

Mi = {l-hio){l + m). (13) 



5 Main results 

In this section we summarize explicit formulae that we were able to obtain 
for the decomposition of the supercritical linear-fractional BGW-processes 
with countably many types. The derivation of these results is given in the 
next section. 



10 



Consider the positive recurrent supercritical case with p > 1 and (3 < oo. 
We demonstrate that the dual reproduction laws are again linear-fractional 

T°° h s- 

/«00 = h l0 + ; , (14) 

with 

^io — — , ^ij — — — , (15) 

% qiP 

1 + m-p A 

m = , = — . 16 

p 1 + m — p 

It turns out that the following remarkably simple formulae hold for the key 
characteristics of the dual branching process 

P = P~\ (17) 

1 oo 

For the Perron-Frobenius eigenvectors we obtain the following expressions 
u = ^"^q- 1 = (q- 1 - 1)(1 + m)(/i - l)" 1 , (19) 

oo 

- = r^E(s H ^- (2°) 



k=0 

We show that the Harris-Sevastyanov transformation results in multivari- 
ate shifted geometric distributions 

m = ; - - , (21) 



where 



1 — q 

hij = _ 3 {hij + mgjiqi - h i0 )) , (22) 

777. 

m = p-l, gj = j—^gj(l-qj). (23) 
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Moreover, we demonstrate that 



P = P, P = jr-Y (24) 



and 



oo 

u = l, v = m^p- 1 - fc (g{H + mp- 1 Hq t g} fc )(l - q). (25) 

k=0 

Theorem 5.1 Consider a linear-fractional BGW-process characterized by a 
triplet (H,g,m). Assume it is supercritical and positively recurrent over the 
state space, that is p > 1 and (3 < oo. Its dual BGW-process and its skele- 
ton are also linear-fractional BGW-processes with the transformed parameter 
triplets (H, g, m) and (H,g,m) with components given by (15), (16), (22), 
@. 

The joint offspring generating function for a skeleton particle of type i 
has the form 

oo 7 \--\oo 7 , 

^MHEtT ^Pl ^rA^+ ^T 1 T" ), (26) 

1 + m — mgs* — (m — mjgt 1 V 1 + m — mgt 1 / 

where 

hij "i(jj(i,h,i,- 



u _ h 



Similarly to the single-type case, we can distinguish in ( 26 ) three components 
but now with dependence: 

• a "reborn" skeleton particle of type i may change its type to j with 
probability hij, 

• independent of % and j a multivariate geometric number of offspring of 
both subtypes, 

• a linear-fractional number of doomed offspring with the fate of the first 
offspring being dependent on ■ 
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The total number of offspring of a skeleton particle of type i has generating 
function fi(s) = Fi(sl, si) of the next form 



/<(*) 



1 + m — ms 



1 — OLi + 



OtjS 



1 + rh — rhs 



where aj = {z^ili—hio) must belong to the interval (0, 1). The corresponding 
mean offspring number is larger than that given by (13): 



Mi = 1 + m + aAl + m) = Mi + (l + m)( h i0 + 



(p- !)(?* - h 



,0, 



6 Proof of Theorem 15.1 

In this section we derive the formulae stated in Section |5l 



Proof of (14). From 



4 /i(sq) _ E.^i Mi?*/ Qi 



Qi 



qi ' l+m-mY™ =1 gjSjqj 



it is straightforward to obtain (14) with (15) and (16). We have to verify 
that gl* = 1 and 

oo 

hio = 1 — hij- 
i=i 



The first requirement follows from (12). The second is obtained from 



Hq* = piHl" - 1* + q*) 
which is proved next. We have (relation (6) in [7]) 



(27) 



H = M 



ni 



1 + m 



Ml*g 



and therefore Ml* = (l + mjHl*, which is (13). Using the last two equalities 



and (11) we find first 
Hq* = 



' Ml* - p(l* - q*) + £ -Ml* 



1 + m 



1 + m 



13 



and then obtain (27). 

Proof of (17). In view of equation ([7]) determining the Perron- Frobenius 
eigenvalue for a linear-fractional BGW-process, to show (17) it is enough to 
verify that 



m 



J^gH™!* = 1. 



n=l 



Observe that according to (15) 

gH fc = ( 

It follows, 



m 



1 + m — p)p 



ir(gH*)q- 



gH n l* 



so that we have to check that 



m 



gH n q* 



(1 + m — p)p n 



n=l 



Turning to (27) we find 



HV = p(H n l t - H™- 1 !* + H"- 1 ^) 



yielding 



(28) 



(29) 



(30) 



(31) 



(p-l^HV^^-q*). 



This and (12) entail (30). 



71=1 



(32) 



Proof of (18). Starting from a counterpart of m\ we find using (29) 



P n=l 



71=1 
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Rewrite (31) as 



nHV = pjntFl 1 - (n - ^H 11 - 1 !' + {n- l)H n_1 q t ^ + /t>H n - 1 (q t - 1*) 
to obtain 



n=l 



n=0 



Thus 



ni 



^gH^-q 



n=0 



" P-l' 



and a counterpart of rt9l) 



Proof of (19) and (20). From (15) we derive H™!* = (H n q*)q • This 



u = (l + TOj/H^p^ETl* 



n=l 



in view of (32) brings (19) 



u 



(l + m)(p- 1) 
p(p- 1) 



oo 

(^HV)q- 1 = (q^ - 1)(1 + m)(Ai - l)" 1 . 



re=l 



On the other hand, a counterpart of (10) together with (28) yields 



v = — — r> P fe gH fe = — > gH fc q 

1 + 777 z — ' 1 + 777. * — ' 

fc=0 fc=0 



Proof of d26b. We have 



/ 4 (s(l-q)+tq)-/ i (tq) 

YL] i Mi (i - + EJli M*fc 



E 7 =i hijtjQj 



+ 



1 + m-m J2kLi QkSkQ- -Qkj-m J2T=i 9ktkQk l + m-m J2T=i 9ktkQk 

1 + m - T77 J2T=1 9kSk - (m - TO) ^fcll ^fc*A 

( Ej°=i Mi<fr) EgLi gfcgfc(i - gfc)) 

(l + 777 - 777 Xl^Li - (TO - TO) J]^=i <Mfc)p(l + m ~ m EfcLl 0***) 
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Replacing the last numerator by m Y^=i 9j s j(^ ~ 1j) J2kLi ^iktkQk an d divid- 
ing the whole expression by 1 — we get 



F(s 



<> = £t 



,-i 



ha + 



1 + m — mgs* — (m — rn^gt* I 1 + rh — mgt* 



and the relation (26) follows 



Proof of (21), (24). Putting t = 1 in (26) we arrive at (21) . Notice 



that according to definition (22) and relations (12), (27) we have 



gl l = 1, HI 1 = 1*. 

Since p is the unique positive solution of 

oo 

m^P^pTl* = 1 

71=1 

and gH n l* = 1 we derive 

oo 

(P-1)£/T n = l- 

71=1 

Thus p = p and 

oo oo 

= rh n/T^glTl* = {p-l)J2 n P~ 



P 



71=1 



71=1 



P-1 
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